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, INTRODUCTION 


Wing tip vortices have long been discussed as a 
primary factor in determining the distribution pattern 
and swath widths of materials ejected from agricultural 
aircraft. Some analytic work has been conducted to de- 
termine the effect of span loading on distribution pattern, 
starting with Wilmer Reed in 1954. (NACA TR //1196, 1953). 
Thousands of distribution patterns have been measured by 
researchers and many attempts have been made to modify 
the pattern of wing tip vortices but none of this work has 
provided a basis for altering either agricultural airplane 
design or agricult ural airplane operations. 

A quantitative assessment of the effect' of span loading, 
both in magnitude and pattern, suitable for use in the field 
or by aircraft designers or modifiers appears to be necessary 
to permit decisions to be made with regard to airplane geome- 
tric characteristics or airplane flight operations. The 
advent of programmable calculators and microprocessor com- 
puters makes it possible to perform calculations which, up 
to this time, have been either difficult or impossible 
because of complexity and length. Accordingly the procedure 
for calculating span wise load distribution as described in 
NACA Technical Report 1071 has been programmed on a program- 
mable calculator, the Hewlett Packard HP-97, as well as in 
BASIC Lanugage. The HP-97 and microprocessors which use 
BASIC Language are readily available and these programs will 
make it possible for either an airplane designer or airplane 
operator to study the effects of span loading on either 
design or operation. 

The procedure of computing a span loading as given in 
NACA Technical Report 1071 has been reduced to a series 
of tables and charts so that the user need only to insert 
the geometric characteristics of his wings, together with 
the operating specifications of his airplane, and a span 
loading for his particular wing can be developed. The 
use of these programs, in conjunction with ASAE Paper //AA 
79-001 will enable operators and designers to explore design 
and operating parameters to determine whether changes or 
modifications are possible that would materially or 
economically effect the swath width and distribution pattern. 
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I 


SYMBOLS 


A 


aspect ratio 



b span of the wing measured perpendicular to 

the plane of symmetry, feet 


c 


wing chord, feet 


'av mean wing chord 


, f eet 


local lift coefficient 


local lift 

qc 


lift coefficient 


total lift 
qS 


6 


rate of change of lift coefficient with 
flap deflection, per radian 


rate of change of lift coefficient with 
wing* angle of attack, per radian 


CiC 


C c 
L 


av 


spanwise loading coefficient for unit lift 


coefficient 


2AG 


G 


spanwise loading 
circulation 


coefficient or dimens ion less 
or _ r_^ 


CfC 


2b 


bV 
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SYMBOLS 


spanwise loading coefficient per radian of 


flap deflection 


M 


m 


S 

V 

W 

W/S 

w 

X 


Mach number 

arbitrary number of span stations 

free-stream dynamic pressure, pounds per square 
foot 

wing area, square feet 

free-stream velocity, feet per second 
airplane weight, lb 
Wing loading, Ib/ft^ 

induced velocity, normal to the lifting 
surface positive for downwash, feet per second 

longitudinal coordinate measured from the 
lateral plane through the quarter chord of the 
wing-root chord, feet 

lateral coordinate measured from the wing-root 
perpendicular to the plane of symmetry, feet 

wing angle of attack, radians 
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SYMBOLS 


a 


V 


da 

d6 


section angle of attack at span station V5 
radians 


lift-effectiveness parameter 



compressibility parame 


ter 




circulation, feet squared per second 


angle of deflection of flap, radians 


dimensionless lateral coordinate 


b/2 


P 

k. 


dimensionless flap span on one wing panel, 
measured perpendicular to the plane of 
symmetry, from the wing root outboard for 
inboard flaps, and from the wing tip 


inboard for outboard flaps 


flap span 
b72 


air density slugs/ft 

ratio of section lift-curve slope at span 


station v to 
Mach number 


27T 


3 


both at the same 




sweep angle of the wing quarter-chord line , 
positive for sweepback, degrees 


wing taper ratio 


root chord 
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SYMBOLS 


SUBSCRIPTS 


f 

t 

r 


pertaining to flaps 
wing tip, 
wing root 


av 

1 

V 


average or mean 

c- 

denoting full wing-chord flaps — = = 1 

c 

pertaining to spanwise station 


CONVERSION FACTORS 


1 meter (M) 

1 sq. meter (M^) 

1 meter/sec. (M/sec . ) 

1 meters /sec . ^ (M/sec2) 
1 Kilogram (Kg) 

1 Kg/M^ 

1 Newton/M^ 


3.281 feet 
10.76 sq. feet 

3.281 ft. /sec. 

3.281 ft. /sec. 

2.205 pounds 
.2048 slugs/ft. ^ 

2.089x10"2 Ib./ft.^ 
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COMPUTATION PROCEDURE 


No attempt will be made to summarize the theory 
or analysis which is the base for computing a span 
load distribution as per NACA Technical Report No, 

1071, Rather, a step-by-step procedure will be tab- 
ulated, which procedure has been programmed for mac- 
hine calculation. 

A . The first step is to determine a span loading 
Cl ♦ c 

coefficient, G = — — ^ This is a dimensionless 

factor which is specified per radian angle of attack (or 
flap deflection angle of one radian) . This coefficient 

is also written as — — = __J!_ , from which it can be 

•^1 2b 

seen that: 

c . c 

1 ^ r 

2b bV 

c . c 


2 




The value of G/6j^ varies with span and is a function 

of wing aspect ratio, taper ratio, sweep and flap span. 

It can be computed from the simultaneous equations (4)‘ 
of Technical Report No. 1071 but it has already been com- 
puted and plotted in Figure 4 of T.R. 1071 for a range of 
aspect ratios, sweep, taper ratios, and flap span. Two 
of these curves. Figures 4(c) and 4(d), for sweep angles 
of 0°, have been used to read values of G/6^ for aspect 

ratios of 6, 8, and 10; taper ratios of .667 and 1.0; and 
flap spans of ,195b, ,556b, and 1.0b. These graphs appear 
as Figure 1 and Figure 2 of this report. G/6j^ values were 

also read for outboard flap spans of .444b and ,805b which 
are the complements of .195 and .556. These values are 
tabulated in Table I through Table V. 
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The first step in the computational procedure is 
therefore to read from Table I through Table V the 
appropriate values of G/6^ for the selected wing and 

flap span at each of the span stations. 

B . The wing specifications of span, taper ratio, 
wing area, and section lift curve slope must be defined. 

^ i ■ 

C . The wing or airplane operating conditions of wing 
loading, air density and speed, in proper units, must 

be stated. 

D . Agricultural aircraft all operate at maximum 
airspeed less than 300 ft./ sec. No correction for 
Mach Number is therefore required. 

E . The wing chord at the spanwise stations corres- 
ponding to values of G/6. must be computed. For a given 
taper ratio , X , the expressions for root and tip chord 
are : 

2 • S 

c = 

b (1+X) 


c = X • c 
t r 


The wing chord at any spanwise station. 


b/2 


is : 


c = c 
y r 


- [w: 


( i-x ) 


F. Since G/6, = 


^1 • ^ 
2b 


. c = 2b (G/6^) 


. c is computed and the wing lift coefficient 

at station y is : Ci . c 

= 1 — I 

1^0 
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G » If a wing operating lift coefficient is chosen, 
the flight speed is: 


V = 


W/S 

(p/2) V‘ 


or 

C 


W/S 


^ (p/2) v2 


H. Thie circulation, F , at each span station is 


<^1„ (P/2) S„ . V2 '^ly • Sy • V 


r = 
y 


pv 


and 


r = 
y 



c 

• y 



2 


V 


A pi 

Ot of Fy 

as 

a function of 

AP/Ay to 

be secur 

ed. 

In order to 

of shed 

vorticity 

wi 

th another win 

must be 

computed 

at 

the same wing 


span will enable 

compare this value 
g, both values of 

C . 

L 


I « The wing lift 
computed as follows: 


coefficient for a 


1 radian 



is 


The wing lift coefficient at any angle of attack is: 

Cli 


"l = 


-• a 


57.3 


The section lift coefficient at any spanwise station 
for a given operating lift coefficient, C is: 

>op 
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J . Wing lift curve slope is: 

dC Ct , 

L = ^1 

da 57.3 

K . The local lift coefficient, at any span 

station ,y, is given by: 
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H . P . 9 7 Program • 


The Hewlett-Packard H.P. 97 programmable calculator 
is a small desk size unit which is programmed in its own 
machine language. 224 program steps are available with 
26 storage registers. Programs are stored on small mag- 
netic cards . 

The Span loading analysis program is listed on pages 
28 and 29 . Data entry instructions and program opera- 

ting instructions are given in appendix A. 

The program is run in two segments using the computa- 
tion forms shown on pages 26 and 27 

The first segment of the program computes the local 
wing chord at station y, ranging from root to tip, at span 
stations as identified from the appropriate Table I through 
V. The wing characteristics are summarized on the first 
computation page, following which 0/6^^ and y/b/2 values 
are listed in the first two columns. The value of Ay is 
the length of each span segment, entered in feet. Succe- 
ssive computations are performed to secure the chord 
length at the midpoint of each span segment. The output 
of the computer is illustrated on page 30 . The computed 

values of are entered in column 5. 

The second segment of the program uses computational 
chart II with the entries as shown on page 31 . The first 

run is an initializing run followed by sequence computation 
of circulation at each span station. A sample computer 
listing with indent if ier co de is shown on page 30 , 
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BASIC PROGRAM 


The span loading analysis procedure has been programmed 
in BASIC language for a Rockwell AIM-65 microprocessor. 

This program should be adaptable to any of the currently 
available microprocessors although it should be carefully 
studied to see whether the syntax of this program is com- 
patible with the microprocessor in question. 

The parameters which describe the particular wing being 
studies, such as gross weight, wing area, lift coefficient 
(corresponding to the desired swath speed), air density, 
taper ratio, aspect ratio, flap span, and so forth, are 
entered manually into the program as per the program state- 
ments. If the value of X in statement 16 is input as zero 
the program computes the swath speed and lists this speed. 
Alternatively a value of 1.0 for X in statement 16 will 
cause the program to compute a flight lift coefficient, 
based upon an input speed. In other words, a choice can 
be made of either selecting an operating lift coefficient 
and the speed will be computed and printed out or the 
swath speed can be input and the operating lift coeffici- 
ent will be computed and printed out. 

A twenty column listing gives sequential values of 
lift coefficient and circulation as a function of span 
position. Final output is the lift coefficient per radian 
for either angle of attack or flap deflection together 
with the lift curve slope of the particular wing. Ad- 
ditional wing geometries can be explored without changing 
airplane parameters by entering the wing geometry para- 
meters in data statements 170, 171, 172, and 173. The 
numbers in these data statements are concurrent values of 
the span station and G/6]^ factor. For the program as ill- 
ustrated in Appendix B the first two numbers in state- 

ment 170 are actual illustrations of the G/6j^ factor of 
.379 at a span station of .05. The numbers as actually 
illustrated in the listed program correspond to values 
for a wing with full span flaps, aspect ratio 8 and a taper 
ratio of .667. These values were taken from Table III. 

In summary, data statements 170 thru 173 are the tabular 
values of the span loading factor as read from table I-V. 

The output of the BASIC program is a function of the 
prxnter bexng used wxth the micropirocessor . A twenty 
column output listing is illustrated on page 34 (Appendix B) . 

Two versions of the BASIC span loading program were 
prepared. These are designated with file names of SPNLA 
and SPNLB respectively. They differ only in the manner 
in which variables are entered into the computer. 
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In SPNLA the airplane and operating parameters are 
written as program statements 17 through 24. This ver- 
sion of the program is useful when a large number of 
wings are to be compared for the same airplane parameters. 

In SPNLB the airplane and operating parameters are 
entered as input statements 17 through 24. With this 
version of the BASIC program every variable must be 
entered each time the program is run but SPNLB gives more 
flexibility in parameter variation. 

The output listing of both SPNLA and SPNLB are the 
same except for listing of input data. A sample output 
listing is shown on page 36. 

WING ANALYSIS 

A large number of wings, with parameters as tabulated 
in Tables VI, VII, and VIII, pages 21 thru 23, have been 
analyzed. Span loading curves for these wings have been 
prepared and are presented in Figures 3&4. These curves, or 
the listing from computer programs, can be used to compare 
span loadings for various wings at various operating con- 
ditions . 

CONCLUDING REMARKS 

Two programs, one for a programmable calculator 
and one written in BASIC for a microprocessor have been 
prepared to compute span loading analysis for wings of 
varying parameters. The span loading as computed by the 
procedures in this report can be used, in conjunction with 
reference 4 to evaluate the effect of span loading and its 
distribution upon the trajectory of particles discharged 
from any point along the wing span. By using the shed 
vorticity, as determined by the method as described in this 
report in conjunction with the dynamics of the discharge 
particles, a prediction of distribution pattern is possible. 
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Table I 

Values of span loading factor G/d^ 
from reference 3 for wings with 
Inboard flaps', ri^ =.195 b/2 


TAPER RATIO=.667 


STATION, %^/2 

ASPECT RATI0=6 

ASPECT RATI0=8 

ASPECT RATIO=10 

.05 

.26 

.235 

.21 

.15 

.225 

.20 

.175 

.25 

.14 

.115 

.09 

.35 ' 

.085 

.068 

.05 

.45 

.06 

.045 

.03 

.55 . 

.035 

.027 

.02 

.65 

.02 

.017 

.015 

.75 

.016 

.015 

.015 

.825 

.013 

.014 

.015 

.875 

.010 

.01 

.01 

.925 

.01 

.008 

.007 

.975 

.01 

.007 

.005 


TAPER RATIO=1,0 


.05 

.245 

.215 

.185 

.15 

.210 

.185 

.16 

.25 

. 125 

.102 

.08 

.35 

• 

o 

00 

.065 

.05 

.45 

.05 

o 

o 

.03 

.55 

.03 

.025 

.02 

.65 

.02 

.017 

.015 

.75 

.015 

.015 

,015 

.825 

.010 

.012 

.015 

.875 

.008 

.009 

.01 

.925 

.008 

. 008 N 

.007 

.975 

.008 

.007 

.005 














































































Table II 

Values of span loading factor 0 / 6 ^ 
from reference 3 for wings with 
inboard flaps^ =.556 b/2 



TAPER RAI 

:iO=.667 


STATION 

ASPECT RATI0=6 

ASPECT RATI0=8 

ASPECT RATIO =10 

.05 

.40 

.343 

.285 

.15 

.395 

.340 

.285 

.25 

.385 

.332 

.280 

.35 

.355 

.302 

.25 

.45 

.30 

.255 

.21 

.55 . 

.235 

.192 

.15 

.65 

.14 

.113 

.085 

.75 

.085 

.065 

.045 

.825 

.055 

.042 

.03 

.875 

.045 

.035 

.025 

.925 

.040 

.030 

.02 

.975 

.025 

.022 

.02 


1.0 0 


TAPER RATIO=1.0 


.05 

.375 

.312 

.25 

.15 

.37 

.315 

.26 

.25 

.36 

.307 

.255 

.35 

.33 

.28 

.23 

.45 

. 29 

.245 

.20 

.55 

.235 

.192 

.15 

.65 

.13 

.107 

.085 

.75 

.06 

.052 

.045 

.825 

.03 

.03 

.03 

.875 

.025 

.025 

.025 

.925 

.025 

.022 

.02 

.975 

.025 

.022 

.02 
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Table III 

Values of span loading factor G/6j^ 
from reference 3 for wings with 
full span flaps Hfsl.O b/2 



TAPER RATI 

>0.667 J 



ASPECT RATI0=6 

ASPECT RAT 10= 8 

ASPECT RATIO=10 

.05 

.445 

.379 

.313 

.15 

.442 

,375 

,307 

.25 

• 

.433 

.365 

.296 

.35 

.419 

.353 

.286 

0 45 

.40 

.337 

.274 

.55- 

.378 

. 319 

.259 

.65 

.347 

.293 

.238 

.75 

.31 

,263 

.216 

.825 

.265 

.229 

.192 

.875 

.225 

,198 

.170 

.925 

.186 

.161 

.135 

♦ 975 

.107 

.09 

.073 


TAPER RATIO=I.O 


.05 

.41 

.345 

,28 

.15 * 

.41 

.345 

.28 

,25 

.408 

.343. 

.278 

.35 

.403 , 

.338 

.273 

.45 

.395 

.33 

.265 

.55 

.388 

.322 

.255 

.65 

.360 

.303 

.245 

.75 

.315 

.273 

.23 

.825 

.275 

.24 

.205 

.875 

.240 

,21 

.18 

.925 

.20 

.173 

.145 

.975 

.12 

.10 

.08 
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Table IV 

Values of span loading factor G/6-i 
from reference 3 for wings with 
outboard flaps, r>^= 0.805 b/2 



TAPER RATIO; 

=0.667 



ASPECT RATI0=6 

ASPECT RATI0=8 

ASPECT RATIO =10 

.05 

.185 

.144 

.103 

.15 

.217 

.175 

.132 

.25 

.293 

.25 

.206 

.35 

.334 

.285 

.236 

.45 

.34 

.292 

.244 

.55 . 

.343 

.292 

.239 

. 65 

.327 

.276 

.223 

.75 

.294 

.248 

.201 

.825 

.252 

. 215 

.177 

.875 

.215 

.188 

.16 

.925 

.176 

.153 

.128 

.975 

.097 

.083 

.068 


TAPER RATIO=1.0 


.05 

.165 

.13 

.095 

.15 

.20 

.16 

.12 

.25 

.283 

.24 

.20 

.35 

.323 

.273 

.223 

.45 

.345 

.29 

.235 

.55 

.358 

.296 

.235 

.65 

.340 

.285 

.23 

.75 

.30 

.258 

.215 

.825 

.265 ; 

.228 

.19 . 

.875 

.232 

.201 

.17 . 

.925 

'.192 . 

.165 

.138 

.975 

. 112 

.094 

.075 
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Table V 



Values of spati 
from reference 
outboard flaps 

TAPER PvAT 

1 loading factor, G/6 
t 3 for wings with 
!, njr=0.444 b/2 

10=0.667 

1 


ASPECT RATI0=6 

ASPECT PvATI0=8 

ASPECT RATIO=10 

.05 

.047 

.039 

.03 

.15 

.046 

.035 

.024 

.25 

.051 

.035 

.02 

.35 

.071 

.053 

.036 

.45 

.099 

.084 

.068 

.55 . 

.152 

.134 

.116 

.65 

.204 

.177 

.151 

.75 

.220 

.193 

.166 

.825 

.203 

.184 

.165 

.875 

.178 

.163 

. 148 

.925 

. 149 

.137 

.125 

.975 

.09 

.076 

.063 

TAPER RATIO=1.0 

.05 

.035 

.031 

.028 

.15 ‘ 

.038 

.03 

.023 

.25 

.05 

.036 

.023 

.35 

.073 

.055 

.038 

1 1 

.45 

.103 

.084 

.065 

.55 

.16 

.136 

.113 

.65 

.218 

.188 

,158 

.75 

,225 

.203 

.18 

.825 

.213 

.196 

.178 

.875 

.193 

,176 

.158 

.925 

.163 

.144 

.125 

.975 

• .103 

.087 

,07 
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TABLE VI 


♦ 


9 


CHARACTERISTICS OF WINGS ANALYZED 


WING 

w/s 

(lb /Ft 2) 

ASPECT 

RATIO 

TAPER 

RATIO 

FLAP 

SPAN 

Of 

WING 

AREA 

(Ft^) 

, WING 

SPAN 
(Ft) 

WEIGHT 

(lbs) .. 

N1 

25 

6,04 

1.0 

None 

326.6 

44,417 



N2 

25 

8.0 

1,0 

None 

326.6 

51.116 

8160 

■N3 - , 

25 

10 

1,0 

None 

326,6 

57,149 

8160 

N4 

25 

6 

.667 

None 

240 

44* 

6000 

■'N5 

25 

8 

.667 

None 

240 

44' 

6000 

N6 " 

25 

10 

. 667 

None 

240 

44’ 

6000 

■■ N7 

25 

6 

.5 

None 

326.6 

44.4' 

8160 

N8 

25 

8 

.5 

None 

326,6 

51.116' 

8160 

N9 

25 

10 

.5 

None 

326,6 


8160 

NIO 

25 

8 

1.0 

None 

450 

60 

11250 

Nil 

25 

7.5 

.667 

None 

450 


11250 

N12 

25 

8 

.667 

None 

450 

60 

11250 

N13 

25 

8 

1,0 

* Inb=,556 

326.6 

mssi 

8160 

N14 

15 

8 

1.0 

Inbn556 

326,6 

51.116 

8160 

N15 

25 

8 

.667 

Inb=.556 

326.6 

IMHHil 

8160 

N16 

15 

8 

.667 

Inb=,55B 

326.6 

51.116 

8160 

N17 

25 

8 

.5 



51.116 

8160 


* inboard of station 






























































TABLE VII 

CHARACTERISTICS OF WINGS ANALYZED 






15 


25 


25 


25 


25 


25 


25 


25 


25 


25 


25 


25 


25 


25 


25 


25 


25 





,667 


,667 


.667 


1.0 
. 667 
1.0 


.667 


.667 


.667 


1.0 


1.0 


0.667 


0.667 


1.0 


1.0 


Inb=.195 

326.6 

Inb=, 195 

450 

Inb=.195 

450 

None 

326,6 

None 

326.6 

None 

450 

1 

* 0. B. Oi. 444 

mam 

O.B. 0.805 

IHlHi 

O.B. 0.444 

326.6 

O.B. 0.805 

326.6 

O.B. 0.444 

450 

O.B. 0.805 

450 


WEIGHT 

(lbs) 




outboard of station 




































































































TABLE VIII 

CHARACTERISTICS OF WINGS ANALYZED 


WING 



ASPECT 

RATIO 


TAPER 

RATIO 


N35 

25 

8 

0.667 

0,B. 0.444 

N36 

25 

8 

0.667 

. ... 

O.B. 0,805 

N37 

30 

8 

.667 

None 

N38 

35 

8 

. ,667 

None 

N39 

40 

8 

.667 

None 

N40 

30 

8 

,667 

Inb=,556 

N41 

35 

8 

,667 

Inb=. 556 

N42 

40 

8 

. 667 

Inb=.556 

N43 

30 

6 

.667 

None 

N44 

35 

6 

,667 

None 

N45 

40 

6 

,667 

None 

N46 

30 

6 

. ,667 

Inb=. 556 

N47 

35 

6 

.667 

Inb=.556 

N48 

40 

6 

.667 

Inb=, 556 







































































































User I nvSl ructions 




INSTRUCTIONS 


Enter Data 


G/61 


Span Station , % °/ 


P ^ S 


Wing Span, b 


Wing Area , S 


Taper Ratio, X 


Gross Weight W 


Air Density, p 
Airplane Cl 


INPUT 

DATA/UNITS 


OUTPUT 

DATA/UNITS 



Airplane Speed, ^^'/sec, input for 



Press A for first run 


Enter data for 2nd run. 


G /&1 


yn 


Ay 


Press B 


Repeat Step III. On last step 

store 1 in Register I, Press B 

Program prints out Cj as last item. 

STO firs’t value of Ay in 

Ay/b/2 


Press C. Using same entries 


ress D. 


Repeat VI for as many steps as nec- 


essary. Last step should be y=l. 0 


& X c„ = 0, last print-out will 







































User Instructions 


APPENDIX A (continued) 
Hewelett Packard 97-7, 
SPAN LOADING ANALYSIS 


INSTRUCTIONS 



INPUT 

DATA/UNITS 


OUTPUT 

DATA/UNITS 


G/5i (From Tables I-V) 


y/ b/2 



Wing Span, Ft. 

Wing Area, Ft.^ 


Taper Ratio, ^^/Cr, A. 


Weight, # lb 





































APPENDIX A (continued ) 


Hewlett Packard 
H. P. 97 


SPAN LOADING ANALYSIS 


Date 

Sheet 


I. Wing, Airplane and Operating Data 


Identification 


Wing Span 

Wing Area 

Taper Ratio ^ 
Operating Weight 

Air Density 

Wing C 

Swath _ 

Root Chord _ 

Tip Chord ^ 
Aspect Ratio 
Flap : 

Type_^ 

Chord 

dC^/da 



ft. 

sq.ft. 


sluRs/cu.ft . 

■ (F6FTT 

—EpE ^ 

ft. 

ft. 


Span 

, dct / dd 


m 

sq.m 

Ic 

k g/m- 
'Constant ) 

Kphr 

m 

m 


% 

,Up to 6^ 


Register 






































































APPENDIX A (Continued) 

PROGRAM LISTING H.P. 97-7 

Span Loading Analysis 


Step 

Key 

Step 

Key 


Entry 


Entry 

001 

LBLA 

055 

0 

002 

DSP 3 

056 

STOI 

003 

SPG 

057 

LBLB 

004 

SPG 

058 

RGLO 

005 

SPG 

059 

PRTX 

006 

P-S 

060 

P-S 

007 

RGLO 

061 

RGLO 

008 

PRTX 

062 

X 

009 

RGLl 

063 

2 

010 

PRTX 

064 

X 

Oil 

RGL2 

065 

PRTX 

012 

PRTX 

066 

ST08 

013 

RGL3 

067 

RGL2 

014 

PRTX 

068 

GHS 

015 

RGL4 

069 

1 

016 

PRTX 

070 

+ 

017 

RGL5 

071 

P=?S 

018 

PRTX 

072 

RGLl 

019 

RGL2 

073 

X 

020 

1 

074 

GHS 

021 

+ 

075 

1 

022 

RGLO 

076 

+ 

023 

X 

077 

P=:=S 

024 

1/X 

078 

RGL6 

025 

RGLl 

079 

X 

026 

X 

080 

PRTX 

027 

2 

081 

STOA 

028 

X 

082 

1/X 

029 

ST06 

083 

RGL8 

030 

PRTX 

084 

X 

031 

RGL2 

085 

P=sS 

032 

X 

086 

STOB 

033 

PRTX 

087 

PRTX 

034 

ST07 

088 

RGLl 

035 

RGL3 

089 

P=sS 

036 

RGLl 

090 

RGLO 

037 

* 

091 

X 

038 

RCL4 

092 

2 

039 

-f- 

093 

tf 

• 

040 

RGL5 

094 

P-S 

041 

* 

095 

PRTX 

042 

2 

096 

STOG 

043 

X 

097 

RGL2 

044 

Vk 

098 

P=^S 

045 

PRTX 

099 

RGLO 

046 

ST09 

100 

X 

047 

P=:=S 

101 

2 

048 

SPG 

102 

-J- 

049 

SPG 

103 

RGLA 

050 

RGLl 

104 

X 

051 

PRTX 

105 

PRTX 

052 

SPG 

106 

STOD 

053 

0 

107 

RGL8 

054 

STOE 

108 

RGLO 
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APPENDIX A (Continued) 


PROGRAM LISTING 
Span Loading Analysis 


Step 

Key 

Entry 

109 

X 

110 

2 

111 

# 

112 

P*S 

113 

RCL2 

114 

X 

115 

RCLE 

116 

+ 

117 

STOE 

118 

PRTX 

119 

RCLI 

120 

X=0? 

121 

GTOE 

122 

RCLE 

123 

v-^s 

124 

RCLI 

125 

-s- ■ 

126 

2 

127 

X 

128 

SPC 

129 

PRTX 

130 

P-S 

131 

ST03 

132 

LBLE 

133 

SPC 

134 

RTN 

135 

LBLC 

136 

RCL3 

137 

5 

138 

7 

139 


140 

3 

141 

-i- 

142 

ST04 

143 

PRTX 

144 

1 

145 

5 

146 

STOI 

147 

RCLi 

148 

RCL3 

149 

-T- 

150 

ST05 

151 

PRTX 

152 

LBLD 

153 

P=?S 

154 

RCL8 

155 

RCL9 

156 

X 

157 

2 

158 

“ 

159 

Ps=S 

160 

RCL5 

161 

X 

162 

ST06 


H.P. 


Step 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 

180 
181 
182 


97-7 


Key 

Entry 

PRTX 

ST08 

RCL7 

RCL8 

PRTX 

RCL9 

+ 

ST09 

RCL8 

ST07 

SPG 

SPC 

RCLI 

X<0? 

GTOl 

R/S 

LBLl 

RTN 

R/S 
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APPENDIX A (Continued ) 


a^G 


.3/3 ST» 


.050 Slot 


. 100 ST02 


Prs 

57, 149 STO0 

IMPUT 

326.600 STOI 

_ DATA 

.667 STD2 


8160. 000 S703 


.002376 S7D4. 


.807 5705 

1 

P:lS 


CSBh 



H.P. 97 COMPUTER LISTING CODE 

for wing N-27 (See Table 


VII) 


PART 

I 

.286 

5700 

.350 

STDI 


eSBB 

0.285 

*** 

32.589 

*** 

5.057 

**f 

5.397 

*t* 


J 0.001 
17.309 




STC'O 


.250 STC'l 


J 


392.574. $$* 
.274 STG0 


0.m 


Pu*T/\T t on 4 
**» FtM?. Nt^-T 
*t7 OUT BaAv.R.D 
~ 6.£06 ««* SPanJ STAT\0M 

5.582 
7.144 ’»♦* 

17.36/ 


293. I6C. 


.216 $ro0 
’50 STCl 
£SBB 
*** 

‘ ♦ 


G.2’t 
2d.6?2 
‘ 5. 

4, ;';-2 

Zl.451 

14.699 

7(4.929 


.t9Z 

.izs 

.050 

0.192 
21.945 
4.973 
4.413 
23.574 
7.105 
746'. 283 


S7D0 

Stoi 

ST02 

QSBE 

$** 

it* 

*** 

*** 

*** 


57.149 

*t*WIKI6 SfW4 FT 

.450 

STJ 

.170 

STO0 

326. see 

***b^lh4& A/ZEA.SaFT 

0.274 

t»oP6 

*875 sm 

0.667 

4** A. TAPER RATIO 



6SBB 

6}6d. 000 

**« &AOSS 

31.318 


0.170 


0.002 

P^Air. demsitv 

5. 829 

ttir 

I9.43J 

$$» 

e.807 

Cl. 

5. 373 


4.859 


6.B56 

(V*, Root ChoRT 

I 2. 859 


3. 999 


4.573 

Tip chord 

15.555 


25. vh 


let. 365 

*** Vv FT/ SEC 

482. 053 


6.942 






774.044 


e.050 

SPAM ST7VT/OM 

.259 

.550 

8700 

870) 

.135 STO0 

0.5/3 

— y 

«« %,(ref 1071) 

0.259 

6SBB 

it* 

. 925 
0.135 

5M/J 

0808 

35.775 


29. 503 



6,742 

*** <j 

5. 60! 

tft 

i5. 430 


5.306 


5.286 


4. ’45 


1,429 

FRoiW (t 

15.716 


3.252 

ttt 

19.266 

»** ARBA. Sl5. RT 

15.004 


25. 4,1 1 

**t 

101.226 

*** A area:) 

555.552 

ttt 

5. 779 


.307 

ST00 ^/S, AT — ^ 

.238 8700 

796. 089 

t*t 

.150 STOI M 

.550 STOI 

.073 

8700 

0.307 

.--nc. 

S/s. 

8.238 

eSBB 

*i* 

.975 

,1.000 

S70J 

STOI 

0800 

35.089 


27. 203 


0.073 

8.344 

6.514 


5.372 

it* 

*•** 

53gr 

t** -CpxP 

5.053 



4,286 

aA f t 

18.57-3 


. - 4. 530 


I8.6t3 

*»* arba . 

15.351 

itt 

1 . 802 
27.850 
5.516 


202.493 

*»» X A AREA3 

644, 7.9;. 

tit 

ft* 



*** 


u* 


WiMGi AT 
oC= \ RADIAM 


30 



APPENDIX A (Continued) 
H.P*. 97 COMPUTER LISTING CODE 
PART II 


.J50STDI .7SQ5TOI 

p:s f-' 

35.089 sm 24.688 ST08 

p:s . 

eSBD ^ 

4gl.735 *** P*A.T 324.869 *** 

9.827 m ^ ~ 33.055 


PR EG 


0.073 0 

0.575 I 
0.050 2 

4.948 3 

0.000 < 

0.000 5 

0.000 f 

0.000 7 

0.000 0 

0.000 5 

4.b'30 0 

:.802 P 
27.0f0 'C 
. p 

scp.eie f 

• J.000 I 


pr^ 

ppft 


.250 STOl 
P?S 

33.832 ST08 

p:s 

GSBD 

445.194 *$* 
16.541 $** 


.350 STOl 

p:s 

32.689 STCS 

p:s 

9SB0 

430.154 **♦ 

]5.04J *u 


.450 STOl 
F73 

31.718 STOP 
PCS 
GSBD 

412.113 
18.041 *»* 


.550 STOl 

p:s 

25.803 ST08 
PCS 
GSBD 

389.545 $4$ 
22.568 *** 


.650 STOl 

p:s 

27.203 ST08 

p:s 

GSBD 


.825 STOl 
.050 ST02 
p;s 

21.945 STOB 

P -*#' 

GSBD 

288. 774 $** 
36.095 *<* 


.875 STOl 
PCS 

19.431 ST 08 
PCS 

, GSBD 
255:692 «*t 
33.082 ♦** 


,.525 STOl 
■ ’ PCS 
. 15.430 STC8 
PCS 
tSBP 

203.043 **» 
52.649 nt 


.975 STOl 
PCS 

8.344 ST08 
PCS 
GSBD 

109.798 $$$ 
93.244 $*$ 


1.000 STOl 
PCS 

0.000 sToe 
PCS 
GSBD 

0.000 *** 

705.758 


57.145 0. 


328.800 

1 

0.887 

P 

6166.00,0 

7 , 

0.002 

4 

0.807 

5 

8.856 

t" 

4.577 

- " ^ - ■ ■ ■ 

8.344 

8 

161.365 

Q 

4,t3P 


1 . 802 

8 

27 . 880 

f 

8.818 

D 

808.010 

8 

1 . 000 

7 

.050 

.<0, OF mvoPoiwT 

STOl ^ OF 

.100 

STO- 


35. 775 

STOS Ca..u X C.A 4 . 

p:s <7 y 

GSBC 

0.088 

JMti 

0.183 


.470.782 


-47C. 782 



CS80 

^70*761 

O.QOO 


IWPUT 


357.564 

3I.S8Z 





APPENDIX B 


BASIC program for SPAN LOADING ANALYSIS 
file name="SPNLA" 


LIST 

1 POKE 41993, 32 

10 PRINT "1071 SPAN LOADING ANALYSIS" 

11 DIM DELY (13) , GDEL (13) 

12 M=1 

14 INPUT "FLAP SPAN =";FSP 

16 INPUT "SPEED COMPUTATION=";X 

17 VFPS=139 

18 WT=6000 

20 AREA=327 

22 RHO=0. 002378 

24 CL=0.8 

26 IF X=1 then 30 

28 GOTO 36 

30 VFPS= SQR (WT/ (AREA*RHO/2*CL) ) 

32 PRINT "SPEED=";VFPS 

34 GOTO 40 

36 CL=WT/ (RHO/2*VFPS'2*AREA) 

38 PRINT "LIFT COEFFICIENT=" ; CL 

40 INPUT "SPAN=";SPAN 

42 INPUT"TAPER RATIO=";TR 

44 INPUT" ASPECT RATIO" :ASPR 

45 CR= ( 2 *ARE A) / ( Sp an * ( 1+TR) ) 

46 PRINT"CR=" ; CR 

47 FOR 1= 1 TO 13 

49 READ DELY (I) , GDEL (I) 

50 NEXT I 

51 DELY=DELY (M) 

52 GDEL=GDEL (M) 

53 K=1-DELY 

54 IF K=0 THEN 80 

PT=TP^PP 

64 CY=CR*(1-(DELY*(1-TR))) 

65 YCL=2*SPAN*GDEL/CY 

66 IF DELY >0.75 THEN 69 

67 YINC=0.1 

68 GOTO 75 

69 YINC=0 . 05 

75 SMCL=SMCL+YCL*CY*YINC*SPAN 

76 M=Mfl 

77 GOTO 51 

80 RCL=SMCL/AREA 

82 MULT=CL/RCL 

84 SLC=RCL/57.3 

90 M=1 
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APPENDIX B (Continued) 


Page two 

91 DELY=DELY (M) 

92 GDEL=GDEL (M) 

102 K=1-DELY 

104 IF K=0 THEN 160 

106 CY=CR* (1-DELY*(1-TR))) 

108 YCL=2*SPAN*GDEL/CY 

110 IF DELY>0.75 THEN 115 

112 YINC=0.1 

113 GOTO 120 

115 YINC=0.05 

120 GAM=YCL*CY*VFPS/2 

125 GAM=GAM*MULT 

126 SGAM=SGAM+GAM’^^INC*SPM/2 

130 PRINT"SPAN STATI0N=" ;DELY 

132 PRINT”CY=";CY 

134 PRINT''CLY=” ; YCL 

140 PRINT"GAMMA=”;GAM 

142 PRINT” ” 

143 M=M+1 

144 GOTO 91 

160 PRINT” ” 

162 PRINT”CL/RADIAN=”;RCL 

164 PRINT”LIFT CURVE SLOPE=”;SLC 

165 PRINT”CT=” ; CT 

166 AVGM=SCAM/ (Span/2) 

170 DATA 0.05,0.379,0.15,0.375,0.25,0.365 

171 DATA 0.35,0.353,0.45,0.337,0.55,0.319 

172 DATA 0.65,0.292,0.75,0.263,0.825,0.299 

173 DATA 0.875,0.198,0.925,0.161,0.975,0.09, 1, 0 

176 PRINT”AVG GAMMA=”;AVGM 

200 STOP 


I 
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BASIC SPAN LOADING ANALYSIS 
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APPENDIX B (Continued) 


BASIC program for SPAN LOADING ANALYSIS 
file name="SPNLB" 


This version of the span loading analysis program is 
the same as SPNLA except for the following statements: 


17 INPUT "VFPS="; VFPS 

18 INPUT "WEIGHT="; WT. 

20 INPUT "WING AREA ="; AREA 

22 INPUT "DENSITY="; RHO 

24 INPUT "LIFT COEFF="; CL 

All succeeding statements are the same as SPNLA. 
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APPENDIX B (Continued) 

COMPUTER OUTPUT LISTING 
for BASIC Span Loading Analysis 

file name=”SPNLB'' 


SPAN STATION^ . 
CY= 6. 52598496 
C'LV- 5. 2880241'; 
396: 325t 



SPAN 

STATION= 


CY= e 

7 2694234 


■ CLV= 

5. 201755 


GAf1HA= 375! 15 




19/i' SPAN L2AU1N0 -AN ^ 

SPAN 

station^ 

ALV’^IS 

CV='6 

. 0138418' 

FLAs-' SPA^4=•:■ 1 i 

CLV= 

4. 963839: 

b F E h D C: F' F U T fl T I n N = b 

GhNNh 

= 343. 40; 

■1 ^ . 

1 5.: 

• 


VFPS='? 139 

SPAN 

STATION= 

UEIGHT=- 6nb,-i 

: CY= 5 

75b’26>02" 

WING AREA-- ■ 3::-7 

■ CLV= 

4. 669295.^ 

DENSITV=? . 90237F " ' 

LIFT C0EFF=? .8 

; GAMMA 

= 309. 29;; 

SPEED= 138. 88835” 

; SPAN 

STATION= 

SPAN=? 51. 116 

r* V r: s 

5665741 

TAPER RATIO=? . 667 

!3LV=" 

4. 2056617 

ASPEC I RATI0='* 8 ■ ’ 

: GAMMA 

= 269. 313 

CR= =' 7. 67512199 ' ' 

SPAN STATI0N= m.=; 

SPAN . 

iTATION = 

CV= 7: 54733121 

.:.CV= '5. 

4Jo7y3i> 

CLV= 13372568 

; CLV= ; 

]. 7217765 

GHMNn= 445. 7193:21 ■ 

' GAMMA: 


SPAN smilQ}i= .1-n--- i : 
CV= ■•'. 29174965 ' 

SPAN ; 

T fi T 1 U fi =•■ 

C^V= 5L 25758588 

1 CLV= ; 

L 0991103 

i3lAMf'1A= 441, Slq.lS . ^ ‘ "I’' 

'■GAMMA= 

189. 342 

SPAN STAT'ION= . 25 

SPAN., i 

*thtion= , 

CV- r", 036-16809 ‘ 

, 'G¥= 5.- 

183281F6 

CLV= -5; 3032672A 

.. C4-V= 1 

7 “7 cr .< 7 A •; 

GAMMA-i , 429. 254755 . ' 

t 

GAMMA= 

105. 843f 

SPAN 9TATI0N= . 35 ' - 



CV= 6, 78058653 . - / 

CL/RAD 

Ihn= 

CLY= 32223811 ' ' ■ 

4, 829 

3694 

QmW-:=- 415. 14227 , 

LIFT C 

br. VC 


. 0842 

b2l68€ 


f ; T = 5 



flVG Gfi 

Mr1fi = 


1 — ' 

99417 


j b’Rh,HK 

T V >v- 


O “t: 
^ !-■ 




36 



X 

0 



FIGURE 1. 



L 


.2 .3 


.4 


.5 


.6 


.7 


.8 .9 


Spanwise 


load distribution factor 


<s 1 


per radian . 


Aspect ratio = 10.0 


FIGURE 2 



BOUND CIRcOlATION - FT^/SEC 


9 



SEMI SPAN POSITION - FEET 































BOUND CIRCULATION -FT^/SEC 



SEMI SPAN POSITION - FEET 











BOUND CIRCULATION — FTVSEC 



SEMI SPAN LOCATION-FEET 


































BOUND CIRCULATION - FT^/SEC 



0 4 8 12 16 20 24 28 32 

SEMI SPAN POSITION - FEET 


« t 




y 700 

w 600 

Ll 


U 300 
on 


200 


O 100 

CQ 


Flap Taper Aspect Area Span W/S ^ Cy Swath 

span ratio ratio sq.ft, ft. Ib/ff^ ^ speed 

mph 


0.556 1.0 8.0 326.6 151.121 25 


0.807 no 


0.556 I 0.667 1 8.0 |326.6 51.12 25 0.807 110 



Figure 3e. Effect of taper ratio 
on boxmd circulation of wings with 
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Figure 3o. Effect of wing loading 
and operating lift coefficient on 
bound circulation of a wing with 
flaps. Aspect ratio=8. 
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Figure 3p. Effect of wing loading 
and operating lift coefficient on 

bound circulation. Aspect ratio=6 
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Figure 4a. Shed vorticity along wing 
80 - span as affected by aspect ratio. 
s= 326.6 ft2. 
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Figure 4c. Shed vorticity as affect- 
ed by inboard flap span. A=.667 
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